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The nature of the resistive transition for a current ap- 
plied parallel to the magnetic field in high-Tc materials is in- 
vestigated by numerical simulation on the three dimensional 
Josephson junction array model. It is shown by using finite 
size scaling that for samples with disorder the critical temper- 
ature T p for the c axis resistivity corresponds to a percolation 
phase transition of vortex lines perpendicularly to the applied 
field. The value of T p is higher than the critical temperature 
for j _L H , but decreases with the thickness of the sample and 
with anisotropy. We predict that critical behavior around T p 
should reflect in experimentally accessible quantities, as the 
I — V curves. 



I. INTRODUCTION 

Dissipation in high-Tc superconductors caused by mo- 
tion of vortex lines has become one of the richest fields 
in the phenomenology of these new materials, (l) From 
the point of view of applications this dissipation is one of 
the principal limiting factors when trying to achieve high 
critical currents. From a purely theoretical point of view, 
due to the high critical temperatures and to the layered 
structure, the physics of the vortex structure becomes 
qualitatively different to that in the low-Tc materials. 
The magnetic field-temperature (H — T) phase diagram 
is not well described by a mean field approximation. An 
estimation of the fluctuation-dominated region around 
the mean field phase transitions in the H — T diagram 
shows that this is quite accessible experimentally. 

One of the new phases of the system that has been 
observed in the High Tc's is the so-called liquid phase. 
This phase is characterized by the fact that the vortex 
lines are easily movable by applying an external force 
on them by putting a current /, i.e., the resistance of the 
system, defined as R — lim/^o V/I is different from zero. 
This liquid phase is found for high enough temperatures, 
whereas for T lower than a well defined value T, the 
resistance becomes zero. Tj is the superconducting tran- 
sition temperature in the presence of the magnetic field. 
The nature of the transition at T depends on the disor- 
der present on the system. When the amount of disorder 
is very small, the low temperature phase is an Abrikosov 
solid phase and the transition at Tj is a first order melt- 
ing transition, with the resistivity having a jump at Tj. 
H For higher disorder a vortex glass phase has been pro- 
posed for T < Tj. pM In this case the transition at T 



is second order and the resistivity displays critical scal- 
ing for T close to Tj. In some cases a crossover due to 
thermal depinning of vortices occurs. |^,|6| 

All these features correspond to the case when the 
magnetic field (which we always supposed applied along 
the c axis of the material) and the external current are 
perpendicular to each other. When the current is par- 
allel to the magnetic field a naive argument would say 
that the dissipation is zero (up to the zero field transi- 
tion temperature T c ) because the Lorentz force on the 
vortices is zero. However, it has been found experimen- 
tally that the transition temperature in this case -which 
we will call T p - is lower than T c . This behavior is related 
to the thermal fluctuations of the vortex lines. 

The point of wether T p coincides or is higher than Tj 
has been highly controversial. It has been found experi- 
mentally that for twinned YBaCuO samples T p turns out 
to be higher than Tj. Q In other words, the supercon- 
ducting coherence is lost in two steps when raising the 
temperature: along the ab plane at T and along the c 
axis a T p . In addition, T p seems to approach T when 
increasing the thickness of the sample, [g Theoretically, 
there have been much discussion about the possibility of 
the loose of coherence in two successive steps, although 
the results are not conclusive. d]l(J Numerical simula- 
tions give support to the picture of one [[Tl-13 or two 
]l4j transitions depending on the model and the param- 
eters used to describe the problem. 

It has been realized that a non-zero dissipation for an 
arbitrary small applied current parallel to the field -i.e., 
a non-zero resistivity- would require the existence of ar- 
bitrary large loops parallel to the ab planes. |I(]Jl4| Ar- 
bitrary large isolated loops have infinite energy, and thus 
the thermal energy which is necessary to create them is 
infinite. However, the superposition of a large amount 
of small loops could generate an arbitrarily large one at 
finite temperatures. The question arises of wether the 
generation of these large loops is a thermal crossover or 
a real phase transition occurring at temperature T p . 

In this paper we will try to clarify these points, by 
performing numerical simulations on the three dimen- 
sional (3D) Josephson junction array (JJA) model, which 
has been used previously to study thermodynamical as 
well as transport properties of high Tc superconductors. 
@,H As the whole problem is too broad, we will con- 
centrate here on the properties of samples with disorder 
and not too high anisotropies, leaving the study of the 
dependence on the amount of disorder for another work. 
pq| It will be shown than in this case -in coincidence 
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with experiments- the 3D JJA model for finite values of 
the thickness has two different transition temperatures Tj 
and T p for the onset of the resistivity in the ab plane and 
along the c axis, respectively. The transition at Ti is -in 
the case of our model and for the values of paramenters 
used- a thermal depinning of vortex lines, but in a real 
materials corresponds to a second order phase transition 
-the vortex glass transition, which has been extensively 
studied before. The temperature T p corresponds to the 
onset of the resistivity along the c axis. We will show |17|| 
that T p is the threshold value of a percolation phase tran- 
sition perpendicularly to the applied held. Moreover, T p 
is a well defined thermodynamical temperature for any 
value of the thickness of the sample, as long as the ab 
plane can be considered as infinite. The precise value of 
T p decreases when increasing the thickness, suggesting 
that the two transitions collapse onto a single one for a 
bulk 3D sample. 

In order to show these properties of the transition at 
T p , we perform finite size scaling on the three 3D JJA 
model. We find that the system exhibits critical behav- 
ior around T„, which reflects, in particular, in a critical 
scaling of the I — V curves. The values of the critical 
exponents are different to those find for the vortex glass 
transition at Tj and agree with those found for the c axis 
V — I curves of real materials. 

The remainder of the paper is organized as follows. In 
the next section we describe the 3D JJA model that we 
use for a detailed investigation of the finite size scaling 
of the transition at T p . In section III a simplified model 
for the percolation transition is introduced. This model 
allows to derive quantitative expressions for the finite 
size scaling that can be checked against the numerical 
simulations, and we do that in section IV, showing that 
the numerical data scale as expected for a percolation 
transition. In section V we focus on the consequences of 
a percolation transition on the experiments, and indicate 
that this consequences are indeed observed. Finally, in 
section VI we summarize and conclude. 



II. DESCRIPTION OF THE MODEL 

For the numerical simulations, we have worked on the 
isotropic three dimensional-resistively shunted-Josephson 



each junction is given by: 
i = I c sin A(f> - 
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junction array model |18| on a cubic lattice. This model 
has been extensively used to derive sensible results for 
thermodynamical as well as transport properties of the 
vortex structure in high-Tc superconductors. When look- 
ing at equilibrium properties it reduces to the 3D uni- 
formly frustrated XY model. 

The 3D- JJA is conformed by a set of Josephson junc- 
tions connecting nearest neighbor nodes of -in our case- a 
cubic lattice. Each junction consists of an ideal Joseph- 
son junction of critical current I c , shunted by a normal 
resistance i?, and a Johnson noise generator, which ac- 
counts for the temperature. Thus, the current i through 



R dt 



(1) 



where A<f> is the difference between the superconducting 
phases on the two nodes connected by the link. A vor- 
tex in a given plaquette of the model is characterized 
by a value of 2ir when following the variation of the su- 
perconducting phase around this plaquette. An external 
current applied between two opposite sides of the sample 
produces a Lorentz force on the vortices of the system, 
and may give rise to dissipation, which is characterized 
by a finite value of the potential difference AV between 
the two contacts. This potential difference is calculated 
in terms of the phases of the system by using the Joseph- 
son equation 



AV 



d (ipx - (p 2 ) 
dt 



(2) 



(pi and </?2 being the phases of the contacts, and (...) 
indicating a time average. 

The system of equations (jl|) are complemented by the 
condition of conservation of currents on each node. The 
equations are numerically integrated in time p9| and dif- 
ferent magnitudes are calculated after thermalization at 
each temperature. 

The boundary conditions (BC) are taken open in the 
ab plane. However, if open BC in the c axis are used, 
there will be a finite force on an isolated vortex at finite 
temperature if the top and bottom ends of the vortex 
are not aligned. The dissipation -that is non-zero even in 
the linear regime- caused by this net force turns out to be 
independent of the thickness of the sample, and in this 
sense, it is only a surface effect. JT5J In order to eliminate 
this spurious surface effect it is crucial to use BC for the 
c direction that assure that each vortex line leaving the 
sample at a given point of the bottom plane re-enters at 
the same point of the top plane. Strict periodic BC on 
the phases ip have this property, however we would ob- 
tain that the voltage difference between top and bottom 
planes is identically zero. 

We use boundary conditions in the c direction that re- 
duce the surfaces effects drastically -although does not 
eliminate them completely. We proceed as follows: we 
numerically integrate the equations of motion for an open 
system one time step, and find the value of the phases in 
the bottom and top plane (pf and (pf . Then the phases in 
this two planes are modified to (pf — ((pf + pf) /2 + (p/2 
and (pf = ((pf + (pf) /2 - (p/2, with (p = (pf - (pf being 
the mean difference between the phases of top and bot- 
tom planes. At this stage the equations are integrated 
another time step. This process guarantees that the vor- 
tex configuration is periodic along the c axis, allowing at 
the same time to calculate the voltage difference when 
we apply a current. The relation between c axis resistiv- 
ity and percolation described in Section IV would not be 
observable if open boundary conditions were used. 
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Flux conservation implies that every flux line going 
into a unit cell of our lattice also goes out of the cell. 
When two vortices go into the same elemental cell we 
cannot tell which one of the two outgoing vortices corre- 
spond to each one of the ongoing vortices. We interpret 
this situation as the meeting of two vortex lines. In a 
real material this corresponds to two vortex lines being 
at a distance lower than the core size of the vortex. At 
high enough temperatures the vortex structure is heavily 
interconnected, and we can follow a vortex line starting 
from one side of the sample and arriving to the opposite 
side. In this case we will say that the vortex structure has 
percolated perpendicularly to the applied field. Due to 
the finite size of the systems used, and to the dynamical 
evolution, percolation is not expected to occur at every 
time, but only at a given fraction of the total time, which 
depends on temperature. We evaluate the probability P 
that there exists a vortex line crossing the system from 
one side to the opposite as a function of temperature. 
This probability P will be one of the most important 
variables in our analysis. 

The ideal Abrikosov (triangular) vortex lattice is frus- 
trated on the subyacent square lattice. The thermody- 
namical properties of the uniformly frustrated XY model 
on a square lattice have strong and non-monotonic vari- 
ations as a function of the field, due to commensurability 
effects. 1 20 This is a spurious effect for us, because we 
are interested in the simulation of systems as close to the 
real samples as possible. In this work we concentrate on 
the case of disordered samples, in which an Abrikosov 
-or any ordered- lattice does not exist at all, even at low 
temperatures. In real samples this may be due to the ex- 
istence of twinned boundaries or point defects, and the 
physics of the low temperature phase may be that of a 
vortex glass or simply a frozen (disordered) set of flux 
lines. 

In the numerical calculations we simulate the disorder 
by considering the 3D JJA system with a random (non- 
correlated) distribution of critical currents of the junc- 
tions. The energy of a vortex in a single loop with one 
Josephson junction of critical current I c is proportional 
to I c , so when taking a random distribution of critical 
currents the links with lowest I c act as pinning centers of 
vortices. If the disorder -i.e, the dispersion of I c - is high 
enough, then the configuration of the ground state will 
be disordered, and not related to the commensurability 
of the vortex lattice parameter and the subyacent square 
lattice. In Fig. [j] we show a top view of the vortex lat- 
tice at zero temperature, as obtained when slowly cooling 
down the system in the presence of an external field of 
value 0.2 in units of flux quanta per plaquette. Two dif- 
ferent results after the annealing process for the same 
distribution of disorder are shown. They correspond to 
a system ofl6xl6x8 sites, with a squared random 
distribution of critical currents ranged between 0.2 and 
1.8, which is the same distribution that we use through- 
out the paper. As can be seen, the vortex structure is 
not the corresponding to an ordered system, and more- 



over the vortex lines are not necessarily straight. Using 
this disordered JJA system we have checked we obtain a 
quite smooth behavior of the quantities we calculate (for 
example the resistivities) as a function of the field. It 
must be keep in mind that the situation for samples with- 
out defects, especially in the case of a triangular lattice 
of Josephson junctions, is different to the one described 
here, [|2l|,P3l and will be addressed in another work. flttl 



III. PERCOLATION TRANSITION AND FINITE 
SIZE SCALING. 

The transition at Ti of our model corresponds to a de- 
pinning of vortex lines from their equilibrium positions 
due to thermal activation. This reflects in the form of 
the resistivity vs temperature curves (current along the 
ab plane) which show a typical thermally activated be- 
havior. jl5|,|| In addition, the process is well described 
by a single vortex model, indicating that collective ef- 
fects are of minor importance. |l5|p3"| 

On the other hand, the transition at T p is a collective 
effect. As we said before, for T > T p there exist some 
paths crossing the sample perpendicularly to the applied 
field, that are free to move when a current parallel to 
the field is applied. We will show now that these paths 
appear at T p due to a percolation phase transition of 
vortex lines. 

We start by introducing a simple model which will al- 
low us to check the proposed percolation phase transi- 
tion by using standard finite size scaling. Let us consider 
the bond percolation problem on a cubic structure of size 
L a b x L ao x L Cl but with one important modification: the 
vertical bonds are supposed to be present with probabil- 
ity one, whereas the horizontal bonds are present with 
probability p. This means that in all cases (for all values 
of p) a rigid squared lattice of lines piercing the sample 
along the c axis (representing the vortex lines at T = 0) 
is present. We will consider the value of p as a function 
of temperature (see below), and look for the probability 
P that there exists a path connecting two opposite sides 
of the sample perpendicularly to the rigid lines. This 
model has a percolation phase transition that can be eas- 
ily characterized in terms of the percolation transition 
of the corresponding two dimensional bond percolation 
problem on the square lattice. For example, if we look 
at the percolation probability Pid on a two-dimensional 
sample of size L a t x L a b, it satisfies (in the limit of large 
L a b) a scaling relation of the form: 



P2D (P2D) = f I (P2D ~ Pc) L a 



-1/u 



(3) 



where / is a universal function. p4f The values of p c and v 
are known to be p c = 0.599 and v = 1.33. The exponent 
v is the one controlling the divergency of a correlation 
length £, i.e., £ ~ {p2D ~ Pc) ■ This correlation length 
measures the typical size of a cluster in the system. 
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The value of P for our model can be read out from the 
solution for the percolation probability of the two dimen- 
sional case, simply noting that two neighbor rigid lines 
are connected if there is at least one horizontal segment 
between them. This corresponds to replace p2D in Eq. 
(||) by 1 — (1 — p) L ". So, we obtain 

P(p) = f(zLli V ) (4) 

with x — 1 — p c — (1 — p) La . Similar scaling expressions 
for other quantities are readily derived. 

We will tentatively apply these scaling expressions to 
the vortex structure. Before showing the results, it is 
important to point out the differences between the two 
situations: the percolation of a real vortex structure is a 
complicated interacting percolation problem, which pro- 
ceeds via thermal fluctuations of vortex lines and ther- 
mally generated loops of different sizes and energy scales. 
Moreover, only paths which conserve the direction of 
the magnetic flux has to be considered. In the simplify 
model, instead, we use a single probability p, which is 
related to the temperature by a Boltzmann factor, i.e., 
p ~ exp (— A/fcT) / (1 + cxp (-A/kT)), with A being a 
energy scale that is roughly given by the energy of a vor- 
tex loop connecting two nearest vortices in the real sys- 
tem. In terms of the temperature, the scaling variable x 
reads 

x = l- Pc -(l + exp(-A/k B T)y Lc . (5) 



IV. NUMERICAL FINITE SIZE SCALING 

We first show in Fig. || the results for the percolation 
probability P obtained for the JJA model as a function 
of temperature -in units of the mean Josephson energy of 
each junction- for an external magnetic field H = 0.2 flux 
quanta per plaquette, which is the value that we used in 
all simulations. The results for different values of L a b 
and L c , as well as the best fitting to equations (^) and 
(||) are shown. The free parameters are p c , A, and v. 
As can be seen in Fig. ||(c), the scaling is quite good in 
spite of the approximations involved, giving confidence 
on the correctness of our ideas. The values obtained for 
the parameters are v = 1.1 ± 0.1, p c = 0.5 ±0.1 and A = 
3.60 ± 0.05. Although these parameters were considered 
as free variables, the relation of p c and v with the values 
corresponding to the 2D model (p c = 0.599, v = 1.33) is 
noteworthy. 

The size effects due to the finite values of L a b and L c 
have very different effects that can be read out from the 
plots in Fig. ^| (a) and (b): when L a {, increases for a given 
value of L c , the percolation transition as a function of 
temperature becomes steeper, indicating that we have a 
well defined percolation transition in the limit L a b — > oo. 
This is the usual finite size effect of a two-dimensional 



second order phase transition. Note that this occurs for 
a fixed (finite) value of the thicknesses L c , and that the 
transition temperature T p -defined as the temperature 
at which there is a jump in the percolation probability 
for L a i, — > oo- is well defined. The threshold value T p 
depends on the thickness L c , as can be seen in Fig. ||(b) 
T p decreases as L c increases. For L C ^S> 1 the form of this 
dependence can be inferred by putting x = in Eq. (|^): 

k B T p ~ A/ln(L c ). (6) 

This indicates a weak decrease of the transition tem- 
perature for thicker samples. This is easy to understand: 
the thicker the sample, the probability of having a vortex 
loop connecting two nearest vortices increases due to the 
fact that there are more places for creating the loop, the 
mathematical expression for this fact being contained in 
equation (||). 

It should be keep in mind that the dependence of the 
transition temperature on thickness is not the usual de- 
pendence of a pseudo-critical temperature on the system 
size, because here we do have a sharp percolation tran- 
sition at any L c (as long as L a b — * oo), and in fact the 
transition is of a two dimensional character. 

From the experimental point of view, it would be inter- 
esting to have a relation between the percolation transi- 
tion and some directly accessible quantity. This quantity 
turns out to be the c axis resistivity of the sample. The 
relation between percolation and c axis resistivity is more 
clearly seen in the model in which the cosine interactions 
of the Josephson junctions are replaced by Villain inter- 
actions. In this case, the model can be exactly mapped 
onto a problem in which -after integrating out a gaussian 
part, related to spin waves- the only degrees of freedom 
that remain correspond to the positions of the vortex 
lines. |2?| Moreover, the voltage between two points A 
and B of the sample is due to the vortex movement and 
can be calculated in the following way: each time a vortex 
line crosses a path joining A and B, the phase difference 
ipAB between points A and B changes in ±2-7r, the sign 
depending on the sign of v x H, where H represents the 
direction of the magnetic field within the vortex, and v is 
its velocity. The mean voltage between A and B is given 
by <fAB/t, where t is the time of measurement. (Nota: 
The argument is asintotically correct only when t — > oo, 
in other case the result may depend on the chosen path 
joining A and B). 

The relation between percolation and c axis dissipation 
can now be understood in the following way: if there are 
mobile vortex lines crossing the sample perpendicularly 
to the c axis, any current along this axis dissipates, i.e., 
percolation of mobile vortex lines is a sufficient condition 
for dissipation. If there are no stationary vortex lines 
percolationg perpendicularly to the c axis, still vortex 
loops can be thermally generated in the sample. These 
loops can cut and reconnect and also blow out and leave 
the sample. The movement of a vortex loop will generate 
a net voltage only in the case in which the loop is blown 
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out of the sample, in any other case only fluctuations in 
the voltage are possible. If a vortex loop grow to infin- 
ity, it implies that it generated during some time interval 
a percolation path across the sample, so percolation is 
a necessary condition for dissipation. Note that the ar- 
gument is not true in the case of using open boundary 
condition along the c direction, because in this case a 
single vortex can give rise to dissipation performing only 
slight displacements from its equilibrium position. Jlf| 

The relation between percolation and c axis resistivity 
was shown before in Ref. JTt| . It was argued that the be- 
havior of p c vs T near the threshold is similar to that of 
S/ (L^ b Ll), where S is the volume of the percolation clus- 
ter (The percolation cluster S is defines as the number 
of elemental segments of vortex lines which are linked to 
(at least) one percolation path) . The argument that led 
to this conclusion was a counting of the vortex lines that 
are involved in the dissipation process, as well as the as- 
sumption of the validity of a viscous motion description 
of the vortex lines. Here we put it in the language of scal- 
ing near the transition temperature. In the same way as 
we derived expression (Q) , we can write a scaled relation 
for the c axis resistivity which reads 

L^p a {T)=g{xL l J b v ) (7) 

with g a new universal function, and t\ and t-i being 
new critical exponents (the value of v in this equation 
should be the same as in Eq. (ji|), if there is only one 
divergent length scale £ ~ (T — T p )~ v close to T p ). The 
c axis resistivity for different values of L a b as well as the 
scaling according to Eq. (0) are shown in Fig. || for a 
system with L z = 8 (the current used to calculate the 
resistivity is always 1/100 of the critical current at zero 
temperature in the corresponding direction). The values 
for the exponents are vi — 0.40 ±0.05 and r 2 = 1.6 ±0.2 
(the value of r 2 is the one used to do the fitting of Fig. 
[|) . The critical behavior of the resistivity in the limit 
L ab — > oo is given by 

P~{T-T P )\ (8) 

where rj is the exponent that determines the asymptotic 
behavior of the function g for large values of its argument, 
i.e., linii^oo g (x) ~ X v . In the limit of very large L ab , the 
dependence on this length should cancel out in Eq. (]?]), 
and we get r\ = vt\ . With the previously found value of v 
and the value of t\ from Fig. ||(b) we find r\ — 0.45 ±0.1. 
The asymptotic behavior of p(T) is shown in Fig. 0(a) 
as a dotted line. 

The close relation between c axis resistivity and per- 
colation can be also seen in the following way. By com- 
bining Eqs. (H|) and ([?]) we obtain 

L2L?p c (T) = g(rHP)) (9) 

In fact, for all the sizes considered, this scaling relation 
is satisfied, as shown in Fig. [|. This figure also shows 



clearly that the c axis resistivity is different from zero 
only if the percolation probability is different from zero. 

The previous scaling theory can be successfully gener- 
alized to account for the possibility of (weak) anisotropy 
in the system. The only change is that now the typical 
energy of an excitation A should be replaced by a value 
which depends on anisotropy. We introduce anisotropy 
in our J JA system by diminishing the c axis mean critical 
current of the junctions and at the same time increasing 
the c axis elemental resistance by the same factor a. In 
a mean field approach, the dependence of A on the pa- 
rameter a is of the form A = Ao / y/a, with Ao the energy 
parameter for the isotropic system. Eq. becomes 

k B T p ~ A A/^ln(L c ), (10) 

making clear that anisotropy gives a much stronger de- 
crease of T p than that obtained when varying L c . 

The validity of Eq. ( f[o| ) is limited at very large L c 
or very high anisotropics, where T p may become close 
to the transition temperature T along the ab direction 
(see below) . For the range of sizes studied, the ab plane 
resistivity when a current of value 0.01 is applied is shown 
m Fig. |. In all cases T ~ 0.6, which is always lower 
than the corresponding value of T p , i.e., we are in the 
region where the previous scaling is expected to be valid. 

V. PHYSICAL CONSEQUENCES OF THE 
PERCOLATION TRANSITION 

The fact that the system has only one relevant length 
scale £ that diverges at T p , has some consequences that 
can be checked experimentally, and in fact, be the key 
point in order to test the adequacy of the theory to the 
experiments. 

From the experimental point of view, the transition at 
T is a vortex glass to liquid transition or a depinning of 
individual vortex lines depending on the strength of the 
disorder. Q In our simulations T always corresponds to a 
thermal depinning. [jl5| The value of T is rather thickness 
independent for the sizes of the sample studied (see Fig. 
j§ and also Ref. JTJ, Fig. 3(b)). 

The dependence of the percolation temperature T p - 
and thus, of the threshold value for the c axis resistivity- 
as a function of the thickness of the sample and the 
anisotropy is given in Eq. (|l0|). As we said above, all 
the results presented here correspond to the case where 
the percolation temperature T p is greater than T, and 
in fact this is a condition for the previous scaling ex- 
pressions to be valid. This is because, in order to have 
dissipation along the c direction we must have a percola- 
tion path across the sample, but it is also necessary that 
this path is free to move. As the path consists of seg- 
ments that go through the horizontal planes, this implies 
that the temperature should be greater than the value 
T for the thermal depinning of vortices in the planes. 
This indicates that percolation is necessary although not 
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sufficient to have dissipation in the c direction. Our pre- 
vious scaling expressions and the relation between p c and 
the percolation probability are valid only if the condition 
T p > Tj is satisfied. 

From the point of view of the numerical simulations, 
it is very difficult to make T p close to Ti by increasing 
the thickness, due to its logarithmic dependence on L c , 
but it is possible by changing the anisotropy a. However, 
when changing the anisotropy, a new feature appears: if 
there is no correlated disorder in the planes, the vortex 
configuration in the limit a — > will be dominated by the 
independent pinning on each planes, with the Josephson 
vortices between planes -which in this limit have an en- 
ergy E — > 0- forming an entangled set of vortex lines. 
If, on the other hand, the disorder is zero, when T — > 
the vortex lines are still straight lines. These different 
possibilities reflect on the relative values of T p and T,-. 

In Fig. [6] we plot the values of Ti and T p as a function of 
the anisotropy for a 24 x 24 x 16 sample . The value of T p is 
given from the onset in the c axis resistivity and from the 
threshold of the percolation transition. The first point to 
be noted is the finite value to which Ti tends for very high 
anisotropics. This value is the depinning temperature of 
a single plane. The two ways of determining T p give 
similar results as long as T p > Tj, and this happens if the 
anisotropy is lower than a critical value a cr . For a > a cr 
the threshold for the percolation transition moves to very 
low temperatures, whereas the onset of the resistivity 
behaves smoothly. The form of the curve T p vs a for 
a < a cr follows the form a -1 / 2 in agreement with Eq. 
©■ 

The situation for a > a cr is not so clear. The fact that 
the vortex loops percolates at very low temperatures indi- 
cates that vortex lines are entangled even for T — > 0, due 
to the disorder within the planes. The dissipation along 
the c direction for T < Ti may be due to the movement of 
vortex lines lying entirely between two neighbor planes. 
In our simulations disorder is a fundamental ingredient 
in order to have T p < Ti. In ordered samples Ti and TL 
collapse onto a single value for high anisotropics, [ p6|j27| 
because the vortices remain straight for T — > 0. In this 
case the threshold of the percolation transition and the 
onset of the c axis resistivity coincide for all anisotropics. 

n 

Experiments performed in twinned YBaCuO samples 
show a decreasing of T p with thickness. This is quali- 
tatively similar to the predictions of Eq. For the 
more anisotropic BiSrCaCuO samples, T p is equal to Ti, 
at least when looking at them by resistivity measure- 
ments, p8| whereas susceptibility measurements suggest 
that T p < Ti, but the value of p c for T p < T < Ti is 
too small to be detected. [^9| This may be similar to the 
behavior in Fig. ^ for a > a cr , although more precise ex- 
perimental as well as numerical work is needed in order 
to make the point clearer. 

In what follows we discuss some experiments that can 
be used to discriminate between a crossover or a real 



phase transition occurring at T p . A percolation transition 
at T p should reflect in the fact that the I—V curves scale 
onto two universal curves for T higher or lower than T p . 
This kind of scaling has been a convincing proof of the 
vortex glass phase transition at T%. 

In our picture of a percolation transition we expect - 
according to scaling arguments near a second order phase 
transition- jl],[| that the I—V curves for different tem- 
peratures plotted as V/I (T - T p ) v vs 1/ (T - T v ) v show 
a universal behavior -i.e., independent of the particu- 
lar value of T-, depending only on wether T is greater 
or lower than T p . The exponents 77 and v can be ex- 
tracted from the results obtained in the previous sec- 
tions. The value at which V/I (T - T p ) v tends for I -> 
when T > Tp should be independent of temperature, so 
V/I ~ (T — Tp) v for low currents. Thus we see that 77 is 
the critical exponent for the resistivity as defined in Eq. 
(||) . We now show that v is the same exponent defined in 
section III. The value of // (T — T p ) u for which the dissi- 
pation starts to be appreciable when T < T p should again 
be temperature independent. This value corresponds to 
the minimum current I cr that is necessary to blow out 
the largest vortex loops perpendicular to the applied cur- 
rent that are present in thermal equilibrium. If these 
largest loops have linear dimensions ~ £, and thus an 
area ~ £ 2 , the energy Ej of this loops in the presence 
of the current I is roughly given by £7 = at; — 6/£ 2 , 
with a and 6 numerical constants, and it has a maximum 
of value Ef = a 2 /AbI for £ = a/267. The energy E 
when there is no current is simply Eq = a£. The loop 
will blow out when Eq > Ef 1 , and we get I cr £ =cte. 
As the typical length £ scales as (T — Tp)~ v we conclude 
that 1 1 (T — T p ) is the appropriate scaling combination, 
thus justifying the use of v for the exponent. It should 
be noticed that the same conclusion can be obtained us- 
ing dimensional analysis, just keeping in mind that 
in our case the relevant dimensionality of our diverging 
clusters is two, instead of three. 

If the proposed mechanism for the transition at T p in 
the 3D JJA model is in fact the one occurring in exper- 
iments, the values obtained for the exponents v and r\ 
in both cases should be similar. Preliminary results in 
YBaCuO samples [[30| indicate that this scaling is rather 
good with values for the exponents that agree with ours. 
Note that the values v ~ 1.1 and 77 ~ 0.45 that we found 
for the exponents are certainly different than those ob- 
tained for the vortex-glass transition along the ab plane 
(v ~ 1.7 and 77 ~ 7). 

In addition to the experimental investigation of the 
scaling of the I—V curves, it would be interesting to 
study the same scaling numerically in the 3D JJA model. 
However, this is a difficult task because of the size ef- 
fects due to the finite value of L a \j. We have done only 
two partial checks of the scaling: we found that the cur- 
rent I cr for which the dissipation starts to be appreciable 
scales as I cr ~ (T p — T) , with the value of v in agree- 
ment with the value obtained independently before . In 
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addition, right at T p and for intermediate values of the 
current, the V — I characteristics show a power law be- 
havior V/I ~ I a with a ~ Tj/u, as expected from the 
scaling analysis, jjj 

The last point we will discuss in order to check our 
percolation theory deals with the behavior of the in-plane 
resistivity p a b as a function of the c axis resistivity p c for 
different values of the external magnetic field. It turns 
out that the form of the p a b vs p c curves for different 
fields can be used to distinguish experimentally between 
a second order phase transition at T p or a crossover due 
to finite thickness of the samples. 

If there were only a unique phase transition at Ti, and 
the effects seen at T p were only due to the finite thickness 
of the sample, the situation would be as follows: The re- 
sistivity would be controlled by an anisotropic correlation 
volume of dimensions £, a b and £ c . Close to Tj we would 
have £afc ~ (T - Ti)~ v and £c ~ (T - T)~ u (notice that 
the numerical factor may be different in both cases), and 
thus at Ti we have ^ a b — * °o and £ c — > oo (we suppose 
that L {, can be considered as infinite in the real sam- 
ples) . The resistivity along the c axis will be zero as long 
as £ c > L c , and due to the very thin samples used, the 
temperature T p at which £ c — L c , turns out to be higher 
than Ti. If the relation £ a f>/£c is field independent, then 
the form of the curves p a b vs p c will not depend on field 
at all. Otherwise, if £ a fc/£c depends on the actual value of 
the field, then p a t vs p c for different fields will not have 
any constraining relation. 

If the system has a percolation phase transition at T p , 
at which a correlation length £ diverges, the situation will 
be the following: the curves p a b vs p c may depend in a 
complicated way on field, except at the point at which p c 
becomes zero. At this point the value of p a b should be in- 
dependent of the field, the coincidence of p a b{p c = 0) for 
different fields can be explained within our proposed per- 
colation transition scheme: the value of the ab plane resis- 
tivity can be calculated using the fluctuation-dissipation 
theorem. For a single vortex this dissipation is propor- 
tional to the typical lateral squared displacement of vor- 
tices divided by the temperature: ~ Ax 2 /T. For the 
whole system the dissipation is proportional to Ax 2 n/T, 
with n the vortex density, which is proportional to the 
magnetic field. So we get p a b ~ Ax 2 H/T. At the tem- 
perature at which p c becomes zero, the displacement Ax 
is roughly given by the distance between vortices, and 
goes as 1 / \/H, and we obtain that p a b(p c — 0) is field 
independent. |3l[| 

Experimental results indicate that the second possibil- 
ity is the one that is realized in high-Tc materials. |3(J It 
is observed that the point at which p c becomes zero cor- 
responds to the same value of p a b irrespective of the field, 
whereas the part of the curves for p c ^ does depend on 
field. This favors our proposal of a real phase transition 
occurring at T p , against the possibility of a crossover due 
to the matching of a correlation length with the thickness 
of the sample. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we presented a description of the c axis- 
resistive transition in the mixed state of high-Tc super- 
conductors for the case of samples where (due to disorder 
and for no so high values of anisotropy) the superconduct- 
ing coherence is lost along the c axis at a temperature 
T p higher than the corresponding to the ab plane. We 
showed that the idea of a percolation phase transition 
perpendicularly to the applied field as the mechanism 
driving the c axis to a linearly dissipative (resistive) state 
is supported by finite size scaling in the 3D JJA model 
as well as transport measurements in YBaCuO samples. 

The transition is a well defined thermodynamic transi- 
tion as long as L a b — * oo, irrespective of the value of the 
thickness L c . However, the threshold value T p decreases 
logarithmically with L c . The analysis of the p a b vs p c 
curves give support to the idea of a real phase transition 
at T p , against the possibility of a crossover due to finite 
thickness with only one single transition at Ti. Also, pre- 
liminary results for the critical scaling of the I —V curves 
are consistent with a phase transition at T p , the critical 
exponents being close to those found in the numerical 
simulations, and certainly different than those found for 
the vortex glass transition. 

As pointed out before, the case corresponding to clean 
samples (that can be numerically simulated using trian- 
gular lattices) gives results qualitatively different to the 
ones described here. In particular, the interplay between 
the transitions for the c axis and for the ab plane gives 
rise, in some cases, to a single first order transition. If dis- 
order is then put into the system, the transition remains 
first order up to a critical value of the disorder. If the dis- 
order is increased further then two separate transitions 
(with T p > Ti if the anisotropy is weak) are recovered. 

The case of very anisotropic samples, where possibly 
T p < Ti, still needs further experimental as well as the- 
oretical clarification. Resistivity measurements in BiS- 
rCaCuO samples have not been able to detect a region 
where p c ^ 0, and at the same time p a b = 0. How- 
ever this could be due to the very high aspect ratio of 
the samples, which prevents an accurate determination 
of p c . On the other hand, our simulations in the 3D JJA 
model suggest that the case T p < Ti is obtained only in 
the case where the vortex lattice -due to the effects of 
disorder- has percolated even for T — > 0. For samples 
without disorder simulations suggest T p = Ti for very 
anisotropic samples. 
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FIG. 1. Top view of the configuration of vortices in a lattice 
of 18x 18x8 for H — 0.2 and a random distribution of disorder 
in the critical currents of the junctions (0.2 < I c < 1.8). Two 
different states after an annealing process are shown. 

FIG. 2. Percolation probability for different values of L ab 
(a) and L c (b) as a function of temperature (in units of the 
mean Josephson energy of the junctions). Panel (c): the 
curves from (a) and (b) scaled according to the given formula 
(see text for explanation) with parameters v = 1.1, p c = 0.5 
and A = 3.6. 

FIG. 3. (a) Scaled c axis resistivity as a function of tem- 
perature for different values of Lab, and (b) the same curves 
plotted against the scaled variable x. The values of the expo- 
nents n and T2 are n = 0.4 and T2 = 1.6. In (a) the limit for 
large L a b -namely p c — (T — T P ) VT1 - is also shown. 



FIG. 4. Scaled c axis resistivity vs percolation probability 
for all sizes quoted in Fig. 2(a) and (b). 

FIG. 5. ab plane resistivity for all system sizes considered 
(as quoted in Fig. 2(a) and (b)). All curves have the critical 
temperature around Tj ~ 0.6. 

FIG. 6. Critical temperatures for the c axis (T p ) and for 
the ab plane (Ti) as a function of anisotropy for a 24 x 24 x 16 
sample. For T p , circles correspond to the onset of the resistiv- 
ity, whereas stars indicate the onset of the percolation. For 
a ^ 10 the onset of percolation moves towards T = 0. 
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